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ABSTRACT 

NGC 3621 is a late-type (Sd) spiral galaxy with an active nucleus, previously detected through mid-infrared 
[Ne V] line emission. Archival Hubble Space Telescope (HST) images reveal that the galaxy contains a bright 
and compact nuclear star cluster. We present a new high-resolution optical spectrum of this nuclear cluster, ob- 
tained with the ESI Spectrograph at the Keck Observatory. The nucleus has a Seyfert 2 emission-line spectrum 
at optical wavelengths, supporting the hypothesis that a black hole is present. The line-of-sight stellar velocity 
dispersion of the cluster is cr* = 43 ± 3 km s"' , one of the largest dispersions measured for any nuclear cluster in 
a late-type spiral galaxy. Combining this measurement with structural parameters measured from archival HST 
images, we carry out dynamical modeling based on the Jeans equation for a spherical star cluster containing 
a central point mass. The maximum black hole mass consistent with the measured stellar velocity dispersion 
is 3 X 10^ Mq. If the black hole mass is small compared with the cluster's stellar mass, then the dynamical 
models imply a total stellar mass of ~ 1 x 10^ Mq, which is consistent with rough estimates of the stellar mass 
based on photometric measurements from HST images. From structural decomposition of 2MASS images, we 
find no clear evidence for a bulge in NGC 3621; the galaxy contains at most a very faint and inconspicuous 
pseudobulge component (M^ > -17.6 mag). NGC 3621 provides one of the best demonstrations that very 
late-type spirals can host both active nuclei and nuclear star clusters, and that low-mass black holes can occur 
in disk galaxies even in the absence of a substantial bulge. 

Subject headings: galaxies: active — galaxies: individual (NGC 3621) — galaxies: kinematics and dynamics 
— galaxies: nuclei — galaxies: spiral 
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1. INTRODUCTION 

In the study of the demographics of supermassive black 
holes, one important unresolved problem is the observational 
census of black holes in very late-type disk galaxies. The 
masses of supermassive black holes in elliptical and early- 
type spiral galaxies appear to be well correlated with th e bulge 
properties of their host galaxies jKo rmendv & Ric hstond 
ll995HFerrarese & MerritdioOOtlGebhardt et al.ii20 00). These 
correlations naturally lead to the question of whether super- 
massive black holes can form in disk galaxies that lack bulges, 
and if so, how the black hole mass in bulgeless galaxies might 
be related to the properties of the host. The nearest exam- 
ple of a bulgeless disk galaxy, M33, does not contain a mas- 
sive black hole, with an extremely tight upper limit of Mbh ^ 
1 500 - 3000 Mq determ ined from s tellar-dynami cal observa - 
tions and modeling (Gebhardt etal.t2001: Merritt etal.l2001h . 
It is not yet known whether the absence of a massive black 
hole is a typical property of bulgeless spirals in general, how- 
ever, because very few late-type disk galaxies are near enough 
for stellar-dynamical measurements to yield such stringent 
constraints on the black hole mass. 

It has long been recognized that many spiral galaxies con- 
tain photometrically and dynamically distinct central star- 
cluster nuclei, with the M33 nucleus being one of th e best- 
studied exampl e s (e.g ., [Walker '1964'; van den Berah 1976'; 
Gallagher et all Il982l |Nieto & Auriere 1982; O'ConneU, 



The Hubble Space Telescope (HST) made it possible 



' Department of Physics and Astronomy, 4129 Frederick Reines Hall, Uni- 
versity of California, Irvine, CA 92697-4575; barth@uci.edu 

^ Department of Astrophysical Sciences, Princeton University, Princeton, 
NJ 08544 

' Hubble Fellow and Princeton-Carnegie Fellow 
The Observatories of the Carnegie Institution of Washington, 813 Santa 
Barbara Street, Pasadena, CA 91 101 



to survey the prope rties of nuclear star clusters in more dis - 
tant galaxies (e.g., iPhillips etal]ll996t ICaroUo etal]ll997h . 
and HST imaging programs have demonstrated that the ma- 
jority of very late-type spirals (Hubble types Scd-Sd) contain 
compact star-cluster nuclei at or very close to their isopho- 
tal centers (Matthews et al. 1999; Bokeretal. 2002). Spec- 
troscopic observations have revealed that the nuclear clus- 
ters in late-type spirals typi cally have stellar v elocity disper- 
sions of - 15 - 35 km s"' (iBoker et al.lll999t IWa lcher et al] 
l2005l) . effective radii of a few parsecs (Boker et al. 2004), 
and dynamical mass es of ~ lO*"- 10^ Mq (Boker et al. 19991 
IWalcher et al.ll2005l) . The spectra also reveal composite stel- 
lar populations indicative of multiple star formation episodes, 
with most clusters containing a population component with an 
age of < 100 Myr (IWalcher et alJl2006l:lRossa et al.ll2006l) . 

The discovery of scaUng relationships between nuclear 
cluster mas ses and host galaxy properties, both in spirals 
iRossa et al. 200 6) and in d warf ellipticals ( Wehner & HarrisI 
2006; Ferrarese et alj2006h . offers an intriguing hint of a con- 
nection between the processes that gov ern the g rowth of nu- 
clear clusters and black holes. Ferrare se et al.l (12006) spec- 
ulate that the formation of nuclear star clusters and massive 
black holes might be mutually exclusive, such that black holes 
are able to form only in galaxies above some critical mass. 
However, the only nuclear star clusters in which the presence 
of a massive black hole can be ruled out at any significant 
level are in the Local Group or within a few Mpc. Opti- 
cal spectroscopic surveys have shown that active galactic nu- 
clei (AGN s) do occur in so me nuclear star clusters in nearby 
galaxies (S eth et aT]l2008ah . so black holes and nuclear clus- 
ters can apparently coexist. Evidence that black holes can oc- 
cur in at least some very late-type disk galaxies comes from 
the detection of a small number of AGNs in Scd and Sd- 
type spirals. The best example is the Sd galaxy NGC 4395, 
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which contains a Seyfert 1 nucleus (iFilippenko & Sargenj 
ll989tlFilippenko & Holl2003T) : it remains the only clear iden- 
tification of a broad-lined AGN in a bulgeless disk galaxy. In 
addition, a few examples of Type 2 AGNs in very late-type 
spirals have been detected re cently in o ptical spectroscopic 
surve ys, such as NGC 1042 dSeth et alJ i2008a: Shields et al., 
l2008h and UGC 6192 (iBarth et al.ll2008l) . In both NGC 4395 
and NGC 1042 the AGNs occur in nuclear star clusters (it is 
not known whether UGC 6192 contains a nuclear cluster). 

A recent Spi tzer spectroscop ic observation of the Sd galaxy 
NGC 3621 by ISatvapal et al.1 ( i 2007) led to the discovery of 
an active nucleus, based on the detection of [Ne V] emis- 
sion lines at 14.3 fim and 24.3 iim. Since photon en- 
ergies greater than 95 eV are required for photoionization 
of Ne"^^ to Ne"^, ordinary H II regions are not expected 
to be significant sources of [Ne V] emission, but a hard 
AGN continuum can easily provide the necessary ionizing 
photons. Emission lines from a range of ionization states 
of neon are observable in the mid-infrared, and the rel- 
ative strengths of these lines are useful as diagnostics of 
the ionization conditions withi n AGN narrow-line regions 
(e.g.,[S pinogHo & Malkan 1 99 21:1 Voitll992l:ISturm et al.l2002l: 
iGroveset aL,2006) . Abel & Satyapall (l2008l) used the results 
of new photoionization models to argue that the strength of 
the [Ne V] emission in NGC 3621 can only be plausibly ex- 
plained by the presence of an AGN, and not by an ordinary 
burst of nuclear star formation. The detection of an AGN in 
NGC 3621 is of significant interest because it is one of the few 
very late-type spirals known to host an active nucleus, making 
it an important target for further observations to constrain its 
black hole mass and AGN energetics. Satyapal et al. (2007) 
note that there is no previously published optical spectrum of 
the nucleus of NGC 3621 suitable for emission-line classifi- 
cation. 

In this paper, we use archival HST images to show that 
NGC 3621 contains a well-defined and compact nuclear star 
cluster. A new optical spectrum of this star cluster is used to 
examine the classification of the active nucleus, and to mea- 
sure the stellar velocity dispersion of the cluster. We describe 
dynamical modeling of the nuclear cluster and the resulting 
constraints on the masses of both the cluster and the central 
black hole. We also examine the structure of NGC 3621 us- 
ing near-infrared images from 2MASS in order to search for 
a bulge component in this late-type galaxy. For the distance 
to NGC 3621, we adopt D = 6.6 Mpc, based on the Cepheid 
measurements of Freedman et al. (2001). At this distance, 1" 
corresponds to 32.0 pc. 

2. IMAGING DATA 

2.1. Archival HST Data 

NGC 3621 has been observed numerous times with HST, 
although many of the observations were of outer fields that do 
not include the nucleus, taken as part of the HST Key P roject 
on the Cepheid distance scale dFreedman et al.l l2001h . We 
found three sets of images that did cover the galaxy nucleus: 
one with WFPC2, one with ACS, and one with NICMOS. The 
parameters for each of these observations are listed in Table 
[U The galaxy was observed with ACS/WFC at two separate 
pointings. One ACS pointing placed the nucleus on the WFCl 
CCD and the other pointing placed it on the WFC2: at both 
positions, identical exposure sequences were obtained in three 
filters. We selected the WFCl pointing since this image in- 
cluded a larger region surrounding the nucleus. The NICMOS 
F190N exposure was taken as the continuum observation for 



TABLE 1 

Hubble Space Telescope Archival Data 



Camera 


Filter 


Exposure Time (s) 


Observation Date 


WFPC2/PC 


F606W 


2x80 


1994-08-17 


ACSAVFCl 


F435W 


3 X 360 


2003-02-03 




F555W 


3 X 360 


2003-02-03 




F814W 


3 X 360 


2003-02-03 


NICMOS/NIC3 


F190N 


6x224 


2007-03-31 



WFPC2/F606W 




NICMOS/F190N 






1" , 


» 













Fig. 1. — //5r images of the nucleus of NGC 3621. Le^ pane/; Aportionof 
the WFPC2/PC F606W image. The rectangle surrounding the nucleus shows 
the position and size of the spectroscopic aperture used in the Keck observa- 
tion, for the l"-wide extraction desciibed in j|3.1l Right panel: A portion of 
the NICMOS/NIC3 F190N image, with a larger field of view displayed. In 
both panels, north is up and east is to the left. 

an F187N (Paschen a) narrow -band image. We do not use the 
F187N emission-line image here. 

The NICMOS and ACS images were retrieved from the 
HST data archives and we use the standard pipeline-processed 
versions of these images. For the WFPC2 data, we use the 
cosmic -ray cleaned and co-added image from the WFPC2 As- 
sociations archive. 

Figure [T] displays the central regions of the WFPC2 and 
NICMOS images. The images show that the galaxy contains 
a very compact and photometrically distinct nuclear star clus- 
ter. Surrounding the cluster is a smooth and nearly featureless 
region of radius ^ 1 ."5 (or ^ 50 pc). At larger radii, dust lanes 
and young star clusters become more prominent in the optical 
images, especially in the ACS F435W band. 

2.1.1. GALFIT Modeling 

To determine the structure of the nuclear clus ter, we use 
the 2 -dimensional modeling package GALFIT dPeng et alj 
120021) . Unfortunately, the nuclear star cluster is saturated in 
the ACSAVFC F555W and F814W images, making it impos- 
sible to use these images to derive structural parameters for 
the cluster The ACS F435 W image is not saturated, but in this 
blue passband the dust lanes and massive stars in the circum- 
nuclear region are more prominent, making this band less eas- 
ily suited to modeling. In the NICMOS image (with a scale of 
0."2 pixel"'), the cluster is unresolved. Therefore, we use the 
WFPC2/PC image for the GALFIT decomposition. WFPC2 
magnitudes listed below are on the Vegamag system, using a 
zeropoint of 22.887 for the F606W filter (Bag gett et aljl997h . 

For the GALFIT modeling, we extracted a 101 x 101 pixeP 
region from the PC image, centered on the nuclear cluster. 
This corresponds roughly to the smooth region surrounding 
the cluster. A larger region would be less suitable for GAL- 
FIT modeling because of the increasing prominence of the 
dust lanes at larger distances from the nucleus. A few patches 
or lanes of dust are visible at the edges of this extracted im- 
age, and we created a mask image to exclude these regions 
from the GALFIT optimization. For the point-spread func- 
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Fig. 2. — Radial profiles of surface brighlness and ellipticity for the nuclear 
cluster in the WFPC2/PC F606W image. "SMA" denotes the semi-major 
axis of the fitted elliptical isophotes. The individual data points represent the 
radial profile measured from the HST image. The long-dashed curve is the 
radial profile of the circularly symmetric GALFIT model, and the solid curve 
is the radial profile of the model fit that was not constrained to have circular 
symmetry. Dotted curves show the radial profiles for the disk and cluster 
separately, for the symmetric model. The middle panel shows the model 
residuals for the radial profile fit (in mag arcsec"^). The slight departures 
from zero ellipticity in the circularly symmetric model are due to random 
errors in the ellipse fitting to the GALFIT model. 

tion (PSF) model, we used the Tiny Tim package ('Krist'1993') 
to create a PSF for the WFPC2/PC camera and F606W filter, 
at 2x oversampling. In GALFIT, the galaxy model is created 
on an oversampled grid, convolved with the PSF model, and 
finally resampled back to the WFPC2/PC plate scale an d con- 
volve d with the WFPC2 CCD charge-diffusion kernel dKristl 
fT993h in order to compare with the observed galaxy image. 

As an initial step, we attempted to model the image as the 
sum of an exponential disk plus a central cluster, using a 
variety of possible models for the cluster based on the vari- 
ous functional forms available in GALFIT (Gaussian, Moffat 
function, Sersic function, and King model). None could ad- 
equately represent the cluster profile; in each case there were 
large systematic residuals. After much experimentation we 
took a different approach, modeling the cluster as a sum of 
several Gaussian components, since a superposition of a suf- 
ficient number of Gaussians can in principle be used to fit an 
arbitrary cluster profile. As a simplifying step, we constrained 
the Gaussian components and exponential disk to be concen- 
tric in the GALFIT model, and found that we were still able 
to fit the cluster profile adequately with this restriction. We 
began with an exponential disk plus a single Gaussian com- 
ponent and added additional Gaussians until the model profile 
ceased to improve significantly with the addition of a further 
component. In our final model we use five Gaussian compo- 
nents plus an exponential disk. 

We created two versions of this GALFIT model. In the first, 
the individual Gaussian components were forced to be con- 



centric and circularly symmetric (i.e., axis ratio = 1). This 
model lends itself most easily to deprojection and dynamical 
modeling, but might not yield an acceptable fit to the cluster 
structure. In the second version, the individual components 
were still set to be concentric, but we allowed each Gaussian 
component to have arbitrary values of axis ratio and position 
angle. 

We used the IRAF^ ellipse task to examine the radial 
profile of the cluster and GALFIT models; this t ask uses the 
elliptical isophote fitting method described by IJedrzeiewskil 
{1987). Figure |2] illustrates the results. We find that the two 
models (symmetric and asymmetric) have nearly identical ra- 
dial profiles. The asymmetric model yields an excellent match 
to the cluster's ellipticity profile (with ellipticity ranging from 

0.015 to 0.38 within the inner r < 2"). Although the sym- 
metric model is obviously unable to match the ellipticity pro- 
file of the cluster, the azimuthally-averaged surface brightness 
profile of the symmetric model is very similar to that of the 
cluster itself, as seen in Figure|2] This makes it possible to use 
a deprojected version of the symmetric model as a simplified, 
spherically-symmetric representation of the cluster structure 
(see ®. 

In the spherical model, the individual Gaussian components 
have FWHMs of 1.8, 4.3, 9.9, 20.9, and 69.7 pixels (with an 
image scale of 0."045 pixel"'). The three narrowest compo- 
nents contribute the majority of the cluster's light. The broad- 
est Gaussian component from the GALFIT decomposition has 
an extent (FWHM=69.7 pixels or 3."1) much larger than the 
visible size of the cluster in the WFPC2 image, and we con- 
clude that it represents light belonging predominantly to the 
galaxy disk rather than to the cluster. As we show below in ^ 
the choice of whether to assign this component to the cluster 
or the disk has a very small effect on the dynamical modeling 
results. The next-broadest component, with FWHM=0"94, 
has a size comparable to the wings of the cluster visible in the 
WFPC2 image (as well as being similar in size to our spectro- 
scopic aperture). We interpret this component as belonging to 
the cluster, although the distinction between cluster and disk 
becomes ambiguous in the outer wings of this component and 
there is no clear way to determine the outermost extent of the 
cluster. 

The F606W magnitudes of these five Gaussian components 
are 17.86, 18.29, 17.97, 18.09, and 17.18 mag, respectively. If 
we represent the cluster by the sum of the first four Gaussian 
components, then its half-light radius is 0."13, corresponding 
to 4.2 pc. Boker et al. (2004) find that 50% of nuclear clus- 
ters in late-type spirals have between 2.4 and 5.0 pc, so 
the size of the NGC 3621 nuclear cluster is very typical for 
this class of objects. The exponential component in the GAL- 
FIT decomposition has a scale length of 23" and a magnitude 
of mpeoew = 9.7 mag, corresponding to V w 10.0 mag based 
on an estimated color of (F606W-y) = -0.3 mag for an Sc 
galaxy template (Kinnev et al.,,1996) . Since the GALFIT fit- 
ting was carried out over only the central 4f'5 x 4 "5 of the 
galaxy, there is no reason to expect this component to be an 
accurate representation of the primary exponential disk of the 
entire galaxy, but its magnitude is surprisingly close to the 
galaxy's total V magnitude of 9.5 mag as listed in the Hyper- 
Leda catalog dPaturel et al.ll2003h . 

' IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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The radial profile plot shown in Figure |2] shows that the 
observed (i.e., PSF-convolved) surface brightness falls by ~ 
5 mag arcsec"^ within the inner r = 1 arcsec of the galaxy 
profile, where the cluster outskirts have essentially merged 
into the surrounding disk. Such a large gradient in surface 
brightness over such a small radius places this object among 
the most solidly detected nucl ear clusters, in com parison with 
the late-type spirals from the iBoker et alj (|2002|) sample. 

2.1.2. Aperture Photometry 

Unfortunately, the saturation in the ACS F555W and 
F814W images precludes their use for GALFIT modeling 
or aperture photometry of the nuclear cluster. If the default 
ACSAVFC GAIN=2 setting had been used for these obser- 
vations, it would still have been possible to perform accu- 
rate aperture photometry on mildly saturated objects. With 
GAIN=2, saturation of the CCD full well results in charge- 
bleeding onto adjacent pixels while still conserving the total 
countrate for the object (Gilliland 2004). However, the NGC 
3621 images were taken with GAIN=1, which does not sam- 
ple the CCD full well capacity, and in this situation the pho- 
tometry of saturated objects is compromised. This leaves the 
ACS F435W, WFPC2 F606W, and NICMOS F190N observa- 
tions as the only usable HST images. The NICMOS image has 
a much coarser plate scale (0."2 pixel"') and a much broader 
PSF than the optical images, and the impact of dust lanes is 
very different between the optical and near-infrared images, 
and as a result it would be very difficult to carry out GAL- 
FIT modeling in a consistent manner for all of these images. 
Therefore, we fall back on a simpler strategy, in which we per- 
form aperture photometry to measure the cluster's brightness 
in each image. 

Under the assumption that the narrowest four Gaussian 
components from the GALFIT model represent the light of 
the cluster, 90% of the cluster's light (in projection) is con- 
tained within a circular aperture of radius 12 pixels, or 0"54. 
We use this as the fiducial aperture radius for the photomet- 
ric measurements, since at larger radii it becomes very dif- 
ficult to distinguish between the outskirts of the cluster and 
the surrounding galaxy disk. For the photometry, we set the 
background annulus close to the cluster, in order to mini- 
mize uncertainties due to the gradient in the fairly flat light 
profile of the host galaxy, and we use a background region 
of 0"6 < r < I'.'Q. For the two optical images, we find 
HST filter magnitudes for the cluster of otf435w = 17.82 and 
inF6()6w = 16.75 mag (both given on the Vega system). Uncer- 
tainties on these measurements are dominated by the choice 
of the inner and outer background radii (rather than photon- 
counting statistics), and we estimate overall uncertainties of 
^0.1 mag on each measurement based on varying the inner 
and outer radii of the background region over wider radial 
ranges. Using the IRAF SYNPHOT package, we can con- 
vert these HST filter magnitudes to Johnson B and V mag- 
nitudes. For a wide range of spectral shapes corresponding to 
SO-Sc galaxy templates ( Kinney etal.ll9 96^. the (F606W-y) 
color index is -0.32 ± 0.02 mag. The ACS F435W filter is 
very close to a Johnson B passband, and SYNPHOT calcu- 
lates a color i ndex o f (F435W-B)= 0.03 ±0.01 mag for the 
iKinney et al.l (Il996h templates. Corre cting for Galactic ex - 
tinction based on E(B-V) = 0.08 mag (ISchlegel et al.ll 19981) . 
we have Mb = -1 1 .65 and My = -12.30 mag for the cluster, 
not including any correction for internal extinction within the 
host galaxy or the cluster itself. 

We carried out photometry using the same aperture and 
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Fig. 3.— Left panel: 2MASS //-band image of NGC 3621. Center panel: 
GALFIT model. Right panel: GALFIT residuals. The image size is 5' X 5'. 
North is up and east is to the left. 

background radii on the NICMOS image, but an additional 
aperture correction was required in order to use photometry 
through an aperture of the same physical size, because the 
NICMOS NIC3 PSF contains significant flux at radii larger 
than 0."54. To estimate the aperture correction, we created a 
GALFIT model of the nuclear cluster at the plate scale of the 
NIC3 camera, using the components from the WFPC2 GAL- 
FIT model, and convolved the model with a NIC3 PSF gen- 
erated by Tiny Tim. We then carried out aperture photometry 
on the simulated cluster image in an 0."54 aperture, and de- 
termined the aperture correction needed to correct that photo- 
metric measurement to one that would correspond to 90% of 
the total cluster light (to be consistent with the optical mea- 
surements). The resulting aperture correction was 0.23 mag. 
Applying the NICMOS photometric zeropoint to the aperture- 
corrected photometric measurement, we find a flux density 
in the F190N filter of U = 2.4 x lO-^*- ergs cm'^ s"' Hz"' at 
19000 A, corresponding to an AB magnitude of 15.43 or a 
Vega magnitude of 13.83 mag. 

2.2. IMASSData 

To determine the overall structure o f the galaxy, we exam - 
ine images from the 2MASS survey dSkrutskie et al.ll2006h . 
Of the three 2MASS bands, the //-band image of NGC 3621 
has the best S/N and image quality, and we used this image 
to create the primary GALFIT model. The image scale is 1" 
pixel"', corresponding to 32 pc pixel"'. While the nuclear 
cluster is unresolved in the 2MASS images, it is still visible 
as a distinct, compact peak in the galaxy center. 

We created a PSF model for the 2MASS image using the 
PSF task in the IRAF DAOPHOT package, combining the 
profiles of 42 individual isolated stars from the image. The 
resulting PSF model is well fit by a Moffat function with 
FWHM = 2"89. 

To construct a GALFIT model of the galaxy, we started with 
a single exponential disk component and added additional 
components as needed to eliminate strong systematic resid- 
uals. A three-component model (convolved with the PSF) 
proved to be the minimal model that yielded an acceptable 
fit, with the three components describing the primary disk, an 
inner disk or possible bulge component, and the nuclear clus- 
ter. We used a Sersic model for the possible bulge component. 
With the Sersic index n left to vary freely in the fit, both n and 
r^ for this component converged to extremely large and unrea- 
sonable values. To better constrain the model, we carried out 
the fit with the Sersic index n of the inner component fixed 
to values 1, 2, 3, and 4, and the model with n = 1 (i.e., an 
exponential profile) yielded the lowest y^. We modeled the 
nuclear cluster as an unresolved source and found an accept- 
able fit with no strong residuals at the location of the nucleus. 
The 2MASS image, best-fitting model, and fit residuals are 
shown in Figure [3] and the galaxy radial profile is shown in 
Figure m Magnitudes for each model component are listed in 
Table 121 For the outer exponential disk, we find a scale length 
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Fig. 4. — Radial surface brightness profile (upper panel) and ellipticity pro- 
file (lower panel) measured from the 2MASS ff-band image using the IRAF 
ellipse task. The individual points show the profile of the 2MASS im- 
age, and the solid curve is the full galaxy model. Individual model profiles 
are shown for the nuclear cluster (dotted line), inner exponential (long-dashed 
line), and outer exponential component (short-dashed line). The nuclear clus- 
ter is modeled as an unresolved source and this component illustrates the ra- 
dial profile of the PSF image. 



TABLE 2 
2MASS GALFIT FITTING RESULTS 



Component 


J 


H 


K 


Primary Disk 


7.61 


6.96 


6.1 i 


Inner Sersic (n = 1) 


12.65 


11.73 


11.47 


Nucleus 


14.14 


13.53 


13.32 



Note. — All results are given in magnitudes, on the 2MASS photomet- 
ric system. 



of h = 45"6 (or L46 kpc) and axis ratio of b/a = 0.48. For the 
inner Sersic component with n = 1, the model yields = 9. "6 
(or 307 pc), and b/a = 0.64. 

Following the same procedures, we created PSF models for 
the 2MASS J and K-hand images, and ran GALFIT on the 
J and K images using the best-fitting //-band model param- 
eters as an initial input to the fit. The fits were very simi- 
lar to the //-band results, and the component magnitudes are 
listed in Table|2] As is typically the case for GALFIT decom- 
positions, the formal fitting errors returned by GALFIT are 
extremely small (< 0.02 mag) and only reflect the statistical 
uncertainties on the best-fit model. The actual uncertainties 
on each component magnitude are dominated by systematic 
errors, typically driven by real structural differences between 
the galaxy's light profile and the idealized GALFIT model 
(such as overall asymmetries, dust lanes, and spiral structure), 
and are probably at least 0. 1 mag. 

The primary exponential disk is the dominant component 
of the galaxy, containing 99% of the total //-band light in the 
galaxy model. Our fitting results demonstrate that the galaxy 
does not contain any substantial bulge component, and pos- 
sibly no bulge at all. At best there is a minor plateau in the 
inner disk surface brightness, corresponding to the inner Ser- 
sic component in the GALFIT fits, that might possibly be as- 
cribed to a faint pseudobulge. The central surface brightness 
of this inner exponential component is 1 mag arcsec"^ fainter 



than the surface brightness of the primary (outer) disk com- 
ponent. 

The inner Sersic component does have some properties in 
common with pseudobulges. It is best fit with an exponen- 
tial profile (n = 1), w hich is typical for pseudobulges in spi- 
ral ga laxies jAndredakis & Sanders 1994; Andredakis et alj 
Il995h . Furthermore, the ratio of its effective radius to the 
scale length of the primary disk is r^/h = 0.21, consistent 
wifli th e mean (r. /h) = 0.22 ± 0.09 found bv lMacArfliur et alj 
(|2003') for a large sample of late-type spirals containing pseu- 
dobulges. With an absolute magnitude of Mh = -17. 4, Mk = 
-17.6 mag, the luminosity of the inner exponential compo- 
nent is within the range spanned b y pseudobulges in nearby 
disk galaxies dCarollo et alj |20011) as well. The inner Ser- 
sic component also has a smaller ellipticity than the primary 
disk, suggestive of a bulge-like feature; this is also seen as the 
gradual increase in ellipticity with radius in the //-band radial 
profile (Figure |4|. Nevertheless, this component is so faint 
relative to the primary disk that it is uncertain whether it rep- 
resents a distinct physical component of the galaxy. Since its 
interpretation is somewhat uncertain, we take its luminosity to 
be an upper limit to the luminosity of a bulge or pseudobulge 
in NGC 3621. Deeper and higher-resolution infrared obser- 
vations would still be desirable, in order to carry out a more 
definitive search for a faint bulge component in this galaxy. 

One possible concern in these model fits is that a small 
bulge, nearly unresolved in the 2MASS images, might be sub- 
sumed within the nuclear cluster component in the model fits. 
However, the small size of this component makes this un- 
likely, because the majority of the light in the unresolved nu- 
clear component is contained within the inner r = 2 pixels (64 
pc). Still, to test for the presence of non-cluster light in the 
nuclear GALFIT component, we compared the 2MASS fit- 
ting results to the small-aperture HST NICMOS F190N mea- 
surement of the nuclear cluster brightness. The NICMOS 
aperture photometry results in = 2.4 x IQ~^^ ergs cm"^ s"' 
Hz~' at 19000 A. Converting the nuclear cluster magnitudes 
from the 2MASS GALFIT models (H = 13.52,/:= 13.32) to 
flux densities, we find f^(l.65 fim)= 4.0 x 10"^'' and /i.(2.2 
/im)= 3.0 X 10"^^ ergs cm"^ s"' Hz"', in each case somewhat 
larger than the NICMOS flux. This suggests that there may be 
some contribution of non-cluster light included in the nuclear 
cluster component in the GALFIT models, which would not 
be surprising given the coarse pixel scale of the 2MASS im- 
ages. Still, the amount of possible "extra" non-cluster light in 
the nuclear cluster model component is well below the bright- 
ness of the inner exponential component in the GALFIT mod- 
els, and the upper limits to the bulge luminosity described 
above are not affected. 

2.3. Stellar Mass Estimates 

The stellar mass of the cluster is a quantity of primary inter- 
est. Ideally, the stellar mass would be determined from HST 
photometry using the widest possible wavelength coverage, 
from the ultraviolet to near-infrared, as well as stellar popu- 
lation fits to spectroscopic data, in order to minimize degen- 
eracies between age, metallicity, and reddening. In addition, 
the stellar population is expected to be composite, containing 
a wide range of stellar ages, as seen in other nuclear clusters. 
The accuracy of our estimates is very limited with only three 
photometric bands available. The lack of photometric cover- 
age blueward of the B band exacerbates the problem, since 
f/-band or near-UV measurements would be needed to ac- 
curately constrain the presence of a young steflar population 
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component in the cluster. Furthermore, the available spectro- 
scopic data (described below) has very low S/N blueward of 
H/3, making it unsuitable for the sort o f detailed s tellar pop- 
ulation modeling that was done by Wal cher et alj ( l2006l) and 
iRossa et aP (l2006l) for other nuclear clusters. 

To estimate the cluster's stellar mass, we used the 
kcorrect package (version 1.4; Blanton & Roweis 2007), 
which performs fits of model spectral energy distributions for 
composite stellar populations with varying amounts of red- 
dening to photometric data in arbitrary filter sets. The code 
us es a basis set of instantaneo us-burst spectra from the models 
of iBruzual & Chariot! (l2003h with the IChabried (l2003h initial 
mass function, and finds the best-fitting match of the spectral 
templates to the photometric data using a nonnegative ma- 
trix factorization technique. To input our photometric data 
points into kcorrect, we correct them for Galactic extinc- 
tion (A = 0.33, 0.22, and 0.04 mag in the F435W, F606W, 
and F190N filters, respectively), and convert to AB magni- 
tudes. The code returns a best-fitting mass-to-light ratio of 
0.47 (in solar units) in the NICMOS F190N band, and a stellar 
mass of M^, = 1.5 x 10^ Mq based on the cluster light within 
the 0"54 photometric aperture. Combining the mass estimate 
from kcorrect with the cluster's V-band luminosity, we 
have a V-band mass-to-light ratio of (M/L)v = 2{Mq /LQy). 

To obtain an additional estimate of the stellar mass from 
the photometric measurements, we applied the results of 
iBell et all (l2003h to derive a K-hand mass-to-light ratio 
based on the optical c olor meas u red fr om the ACS and 
WFPC2 images. The iBell et al.1 (l2003l) M/L values are 
bas ed on th e PEGASE spectrophotometric evolution mod- 
els (iFioc & Rocca-Volmerange 1997), assuming an exponen- 
tially declining star-formation rate. The assumed stellar initial 
mass function (IMF) is a modified Salpeter (1955) IMF with a 
smaller contribution from low-mass stars, resulting in a 30% 
reduction in mass relative to a Salpeter IMF for a given lu- 
minosity. Bell et al. estimate that their M/L ratios may have 
^ 20% random uncertainties and ~ 25% systematic uncer- 
tainties due to dust extinction and differences in galaxy ages 
and star-formation histories. Uncertainties in the IMF are 
likely to contribute at least this much systematic error as well. 

Transforming the optical photometry to Johnson filters, we 
find (B-V) « 0.6 mag for the cluster. For this color index, 
the Bell et al. results predict (M /L)k = 0.75 in solar units. 
Applying this mass-to-light ratio to the point source com- 
ponent in the 2MASS /T-band model gives a stellar mass of 
3 X 10^ Mq. As discussed above, the comparison between 
NICMOS and 2MASS photometry of the nuclear cluster sug- 
gests that the 2MASS K magnitude might be too bright by 
^ 20% due to inclusion of non-cluster light. The Bell et 
al. prescriptions yield a V-band stellar mass-to-light ratio of 
(M/L)v = 1.4(M0/L0,v), and combining this with the abso- 
lute V magnitude determined from the WFPC2 data implies 
= 1.0 X 10^ Mq. From all of these results, our best esti- 
mate of the cluster's stellar mass is (1 - 3) x 10^ Mq. 

These stellar masses are highly uncertain due to the limited 
spectral coverage of the HST data. Even with a more complete 
dataset, however, there are potentially large systematic uncer- 
tainties in st ellar masses obtained by these techniques. As 
discussed by iKannappan & Gawiserl (|2007|) . different meth- 
ods for estimating galaxy stellar masses (i.e., from spectral 
energy distribution model fitting with different evolutionary 
synthesis models, or from the Bell et al. M/L calibration) can 
result in systematic differences of at least a factor of 2, and 
possibly more than a factor of 3, even when the same IMF 



is assumed. Allowing for uncertainty in the IMF, a system- 
atic uncertainty of more than a factor of 4 in stellar mass is 
certainly plausible, particularly in view of the complex star 
formation history of a typical nuclear star cluster. Ultimately, 
the main conclusion that we can draw from these results is 
that within the substantial uncertainties, the stellar mass de- 
rived from the HST photometry appears to be consistent with 
the cluster's dynamical mass (as described below in ® . 

3. SPECTROSCOPIC DATA 
3.1. Observations and Reductions 

A 600 s exposure of the nucleus of NGC 3621 was obtained 
at the Keck-II telescope on the night of 2008 March 2 with 
the E chellette Spectrograph and Imager (ESI; ISheinis et al.l 
l2002h . An 0."75-wide sHt was used in ESI echelle mode, 
and the slit was oriented at the parallactic angle (PA=9°). 
The instrumental dispersion with this slit width is (T,- w 22 
km s"', and the pixel scale in the dispersion direction is 
1 1.5 km s"' pixel"', with total wavelength coverage of 3850- 
11000 A over 10 echelle orders. We do not have a direct 
measurement of the seeing during the exposure (taken at air- 
mass 1 .62), but the seeing in standard star observations on that 
night was ^ 0"8- 1"2. Extracted spectra were wavelength- 
calibrated using observations of HgNe, Xe, and CuAr com- 
parison lamps, and flux-calibrated with an observation of the 
standard star Feige 34. Error spectra were extracted and prop- 
agated through the same sequence of calibrations. 

We performed the spectral extractions twice using differ- 
ent extraction widths. To obtain a spectrum of the nuclear 
cluster, we extracted the spectrum using an extraction width 
of 1", and a background region of 2"5-5."0 on either side of 
the nucleus. This small-aperture spectrum did not yield sat- 
isfactory results for the emission lines, however The emis- 
sion lines are extended over a region of several arcseconds 
along the slit, and the emission-line surface brightness is not 
centrally concentrated. Consequently, the small-aperture ex- 
traction included only a small fraction of the total emission- 
line flux in the slit, and the close-in background regions had 
nearly the same emission-line surface brightness as the ex- 
traction window. In the resulting spectrum, the emission lines 
had erroneously low fluxes due to this background subtraction 
problem. To obtain a better emission-line spectrum, we per- 
formed another extraction using a 5" extraction width, with 
the background regions close to the edges of the slit where 
the emission-line surface brightness is lower We use the nar- 
row extraction in ^3.2l below for measurement of the nuc lear 
cluster's velocity dispersion, and the wider extraction in ^3.31 
to examine the nuclear emission-line ratios. 

3.2. Stellar Kinematics 

The velocity dispersion of the nuclear star cluster wa s mea- 
sured with the direct-fitting routine described by Bart h et aP 
(2002), using spectra of four stars of spectral type G8III- 
K2III observed during the same run. We performed fits sep- 
arately to spectral regions around Mgb (5100-5450 A) and 
the Ca II triplet lines (CaT; 8400-8750 A). The fitting routine 
models the galaxy spectrum as the sum of a stellar spectrum 
broadened by a Gaussian kernel and a power-law featureless 
continuum. The velocity-broadened spectrum is multiplied by 
a quadratic polynomial that can account for both reddening 
and minor differences in flux calibration between the galaxy 
and template star. The resulting spectrum is then fitted to the 
data, with the velocity dispersion, power-law slope and nor- 
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Fig. 5. — Spectral regions used for measuring the stellar velocity disper- 
sion. The red curve is the best-fitting velocity-broadened template star spec- 
trum. 

malization, and multiplicative polynomial parameters allowed 
to vary as free parameters in the fit. The measurements for 
the blue and red spectral regions yielded identical results of 
(T = 43 ± 3 km s~', where the uncertainty represents the sum 
in quadrature of the fitting error on the best-fitting template 
and the standard deviation of the velocity dispersions derived 
from all template stars. Figure |5] shows the template fits to 
the spectrum of the nuclear cluster The Gaussian-broadened 
template spectrum fits the profiles of the strong Ca II lines 
well, indicating that the possible contribution of higher-order 
moments to the line-of-sight velocity profile is relatively in- 
significant for the integrated spectrum of the cluster. 

For the blue spectral region, differences in the the Fe/a 
abundance ratio between the galaxy and the template stars 
can potentially lead to biased results, because the strongest 
features immediately redward of Mgb are Fe features such 
as Fe5270. The direct-fitting routine requires the featureless 
continuum dilution to vary smoothly with wavelength, and it 
cannot properly account for strong differences in line strength 
between adjacent Fe and Mg features (see Barth et al. 2002 for 
further details). However, in this case we found that the fit- 
ting results were not sensitive to whether the Mgb lines were 
included or masked out in the fit, and Figure |5] shows that 
the routine yields a good fit over the entire 5100-5450 A re- 
gion. As an additional check, we also measured the velocity 
dis persion using a n indep endent direct-fitting code described 
by iGreene & Hoi (l2006ab and found very consistent results 
(44 ± 2 and 43 ± 2 km s~' for the Mgb and CaT regions, re- 
spectively). We adopt (T* = 43 ± 3 km s"' as the best estimate 
of the line-of-sight stellar velocity dispersion. 

The ESI slit length of 20" is too short to obtain a pure sky 
spectrum simultaneously with the galaxy observation, and it 
is unavoidable that the background regions used for sky sub- 
traction contain light from the galaxy disk. Consequently, the 
background subtraction employed in the spectral extraction 
must be subtracting off some galaxy light with a different ve- 
locity dispersion and possibly a different mean velocity than 
that of the cluster itself. We tested the magnitude of this effect 
using the data from the ESI echelle order containing Mgb and 
Fe5270. For this order, we added the extracted sky spectrum 
back into the extracted cluster spectrum, to obtain a spectrum 



of the cluster without sky subtraction. We then ran the veloc- 
ity dispersion fitting routine on this spectrum. The result was 
an increase of only 1 km s"' in the measured dispersion, in 
comparison with the value measured from the sky-subtracted 
spectrum. Thus, we conclude that the sky subtraction is not 
adversely affecting the velocity dispersion measurement. In 
the spectral order containing the Ca II triplet lines, the back- 
ground spectrum is dominated by night sky emission lines, 
making it impossible to carry out a similar test for the Ca II 
dispersion. 

Even when background subtraction is applied during the 
spectral extraction, the spectrum includes some residual light 
from the underlying galaxy disk. We used the results of the 
GALFIT decomposition of the WFPC2 image to estimate the 
fractional contribution of galaxy disk light to the ESI spec- 
trum in the V band. We took the GALFIT models for the nu- 
clear cluster and disk separately, convolved each with a Gaus- 
sian kernel of FWHM = 1 "2 to represent the likely seeing dur- 
ing the Keck observation, and measured the flux from each 
component within an aperture corresponding to the ESI slit 
(0"75 X 1"). Within this aperture, the disk contributes 26% of 
the total light. However, the background subtraction used in 
the spectroscopic extraction reduces this disk contamination 
substantially. The background regions used were 2."5-5."0 on 
either side of the nucleus. We used the PC image itself to 
estimate the flux level in the background regions, since the 
GALFIT model does not include real structure such as dust 
lanes that are present in the actual sky background regions. 
We convolved the PC F606W image with the same Gaussian 
kernel and measured the counts in the exact background re- 
gions used in the ESI extraction. After subtraction of the 
mean background level, the remaining contribution of galaxy 
disk light to the ESI spectrum is only 10%. From the results 
of the background-subtraction test described above, we con- 
clude that this small residual contamination of disk light in 
the ESI spectrum should not have a significant effect on the 
measured velocity dispersion of the cluster. 

It is important to note that our measured value of (J^, is 
essentially the second moment of the cluster's integrated 
line-of-sight velocity profile, which may contain a contribu- 
tion from rotation as well a s from random stellar motions 
within the cluster. ISeth et af] (l2008b ) have found from high- 
resolution integral-field spectroscopic data that the nuclear 
cluster in NGC 4244 (an edge-on Scd galaxy) is rotating at 
30 km s"', and if nuclear clusters are formed in situ by in- 
flowing gas in the host galaxy disk, rapid rotation might be an 
important feature of nuclear clusters in general. 

3.3. Emission-Line Spectrum 

To measure the emission-line strengths and ratios, we car- 
ried out starlight subtraction on the wide-aperture spectral ex- 
traction using similar techniques, this time using a linear com- 
bination of a G8III-K3III giant and an AOV star spectrum to 
represent the stellar continuum. The starlight-subtraction fits 
were carried out over the ranges 4600-5300 A and 6000-7000 
A, with the wavelengths of the narrow emission lines masked 
out in the fit. The results are shown in Figure |6] The starlight 
subtraction works reasonably well but results in an imperfect 
fit to the Ha and H/3 absorption lines, leaving some spuri- 
ous residual absorption in the starlight-subtracted spectrum. 
This residual absorption is the main source of uncertainty in 
the emission-line flux measurements described below, partic- 
ularly for the very weak H/3 line. 
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Fig. 6. — Starlight subtraction in the spectral region surrounding the H/3 
and [O III] lines (upper panel) and the Ha and [N II] lines (lower panel). In 
each panel the spectrum of the nuclear cluster is shown in black, the best- 
fitting starlight model in red, and the starlight-subtracted spectrum in blue. 

Emission line fluxes were measured by direct integration 
of the starlight-subtracted spectrum, with the results listed 
in Table [3] For Ha and H/3, we attempted to correct for 
the starlight-subtraction residuals by manually fitting a spline 
to the residual absorption profile and subtracting it from the 
spectrum; we performed several trial subtractions and the un- 
certainty range in the measured fluxes reflects the range of 
values measured in different trials. Although some line fluxes 
are highly uncertain, the results provide a clear emission-line 
classification for the nucleus. Figure [T] shows the location of 
the NGC 3621 nucleus in a line-ratio diagnostic diagram, in 
which it falls in the main Seyfert branch of the diagram. Thus, 
we classify NGC 3621 as a Seyfert 2 galaxy, and the optical 
spectrum adds further evidence that this galaxy contains an 
AGN. 

From inspection of the two-dimensional ESI spectrum, the 
emission-line properties change systematically at a distance 
of roughly 3 "5 from the nucleus. At this radius, the [O IIIj/Ho! 
and [N Ilj/Ha ratios appear to drop substantially. If this 
change in emission-line ratios marks the boundary of the 
AGN-dominated narrow-line region, then the radius of this 
region is ^ 1 10 pc. Assuming a circularly symmetric narrow- 
line region on the plane of the sky, the ESI spectrum in- 
cludes only ~ 10% of the total area of the emission-line re- 
gion. Within this region, the emission-line surface brightness 
appears nearly constant in the ESI spectrum, thus our spec- 
trum probably contains only about 10% of the total narrow- 
line flux of the AGN. It would be very useful to obtain a nu- 
clear spectrum with an integral-field spectrograph, in order to 
better map out this emission-line region and obtain improved 
estimates of the total AGN luminosity. 

The emission-line widths are unresolved in the ESI spec- 
trum. We find FWHIVI w 52 km s"' for arc lamp lines ob- 
served through the same 0"75 slit width. The raw linewidths 
(not corrected for instrumental broadening) are 5 1 km s"' for 
[O III] A5007 and 57 km s"' for [N II] A6583, so the intrinsic 
linewidths must be substantially smaller than the instrumental 
resolution. Such narrow lines are exceptionally unu sual for 
Seyfer t 2 nuclei, even in low-mass host galaxies. Barth et al.l 
( |2008|) used this same ESI setup to measure the linewidths 



TABLE 3 
Emission-Line Measurements 



Line 


Flux 


H/3 


2.6± L3 


[O III] A5007 


22.4 ±0.4 


Ha 


8.9 ±2.4 


[N II] A6583 


13.0± L8 


[S II] A6716 


2.0 ±0.5 


[S II] A6731 


1.8±0.5 



Note. — Emission-line fluxes are given in units of 
10"'* ergs cm"^ s"', and ai'e corrected for Galactic red- 
dening of E(B-V) = 0.08 mag i Schlegel et al. 1991). 
The flux scale has an unknown overall scaling uncer- 
tainty due to slit losses. 




log ([N II] / Ha) 

Fig. 7. — Emission-line diagnostic diagram illustrating the AGN classi- 
fication of the NGC 3621 nucleus. Contours represent SDSS emission-line 
galaxies from Kaufl'mann et al. (2003). The dashed and dotted lines represent 
the classification boundaries of Kewle v et aU t2006) . 

for a sample of low-mass Seyfert 2 galaxies selected from 
the Sloan Digital Sky Survey, and the smallest [O III] width 
found was FWH]VI=66 km s"' in the late-type spiral UGC 
6192. NGC 3621 is so nearby that the ESI spectrum sub- 
tends only a small portion of the narrow-line region, and the 
linewidth might be substantially different if measured through 
a larger aperture. 

The internal reddening toward the narrow-line region can- 
not be determined accurately, because of the very uncertain 
H/3 flux measurement. From the Balmer line fluxes (cor- 
rected for Galactic extinction) we find Ha/H/?= 3.4 ± 1.9. 
This implies a possible range for the internal reddening from 
E(B-V) = up to E{B-V) = 0.6 mag, assuming Case B re- 
combination for the intrinsic Balmer-line ratios. 

4. DYNAMICAL MODELING 

As an initial approximation of the cluster's total mass, we 
use the virial theorem for a spherical, isotropic, and isolated 
system: 

Mvi, = /^. (1) 

Here, / is a dimensionless factor that depends on the density 
profile and velocity dispersion profile of the cluster, and a^, is 
the line-of-sight velocity dispersion averaged over the density 
profile of the cluster. The numerical value of / is typically 
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~ 10 (e.g.. IWalcher et ani2005l) if the effective radius r^, is 
used as the cluster radius parameter. Taking re « 4 pc and 
(7* = 43 km s"', this yields Mvir « 1 .7 x 10^ Mq. 

Though the virial theorem provides a rough estimate of the 
cluster mass, the precise value of the factor / depends on 
the density profile and velocity dispersion profile of the clus- 
ter. Additionally, the value of / would be affected if there 
are multiple components that significantly contribute to the 
mass of the cluster, such as a point mass in the form of a 
central black hole. In this case, the outskirts of the cluster 
blend smoothly into the galaxy disk, and the outer radius of 
the cluster is not well determined. For these reasons, in order 
to derive a more accurate cluster mass, as well as to determine 
constraints on the mass of the black hole, we compute dynam- 
ical models for the cluster based on the Jeans equation. This 
formalism allows us to model the structure of a star cluster 
with a two-component mass distribution including stars and a 
central point mass. 

Similar to the virial theorem method, the total cluster mass 
as derived from the Jeans equation depends on its outer extent. 
However, determining the extent of the light that is gravita- 
tionally bound to the cluster is somewhat complicated by the 
fact that it is embedded in a disk. Particularly near the out- 
ermost radii, it is not entirely straightforward to distinguish 
between the light from stars bound to the cluster and the light 
due to the stars in the disk. In the GALFIT model for the 
cluster, 90% of the total cluster light in the multi-Gaussian 
model is contained within a projected radius of 0."54, or ^ 17 
pc. To conservatively model the mass bound to the cluster, 
we use Tout = 17 pc as the fiducial outer radius of the cluster. 
All results described below from the dynamical modeling re- 
fer to the mass enclosed within this radius. We also note that 
the results described here refer to the mass of the nuclear clus- 
ter; the mass contributed by the galaxy disk component within 
this fiducial radius is not included in the dynamical mass re- 
sults quoted below. Within a projected radius of 17 pc, the 
total light from the four GALFIT components corresponding 
to the cluster is equivalent to L = 7.8 x 10^ l^Qy- 

The Jeans equation for the radial velocity dispersion, CTr, is 



dr 



GM(r) 



(2) 



The anisotropy of the stellar distribution is defined as /3 = 1 - 
ag I aj, the total dynamical mass enclosed within radius r is 
M(r), and the three-dimensional luminosity density is p+(r). 
We use the boundary condition of p^,a^ at the outer radius 
of the cluster, which we take to be 17 pc as described above. 

To determine p*(r), we deproject the light distribution of the 
cluster via an Abel transformation. We use the results of the 
circularly symmetric GALFIT decomposition of the WFPC2 
image as the model for the observed light distribution, assum- 
ing that the four narrowest GALFIT Gaussian components be- 
long to the cluster. These Gaussian components are individ- 
ually deprojected by Abel transformations, and the resulting 
p*(r) is the sum of these four deprojected functions. 

To compare to the measured velocity dispersion of the clus- 
ter we integrate along the line of sight through the cluster: 



a\R): 



■dr. 



(3) 



h{R) Jr \ r^J Vr^ 

where R is the projected radius on the plane of the sky and 
h{R) is the surface brightness at R. We carry out the integra- 
tion up to a radius of ro^t = 17 pc, since this is the assumed 
outer radius of the cluster Given the multi-component light 



distribution determined from GALFIT, equations|2]and|3]can- 
not be solved analytically, so to obtain we first numerically 
integrate the Jeans equation to obtain aj, and then integrate 
this solution along the line of sight using equation|3]to obtain 

a\ 

Finally, to obtain the integrated line-of-sight velocity dis- 
persion for comparison with the observed value, we average 
the line-of-sight velocity dispersion over the cluster, weighted 
by surface brightness: 

2tt J a^iR)h{R)RdR 



(4) 



Itt J UR)RdR 

In the form above, the Jeans equation assumes spherical 
symmetry; given the small observed ellipticity of the light dis- 
tribution of the cluster, as well as the fact that the circularly 
symmetric GALFIT model fits the cluster's radial profile well 
overall, this provides a good approximation to the dynamical 
mass. Further, this approach assumes that the cluster is non- 
rotating. Though the data do not give us any information on 
the rotation of this cluster, the small ellipticity can be inter- 
preted as a small amount of flattening of the cluster due to 
rotation. 

We model the mass distribution, M(r), in terms of a compo- 
nent proportional to the stellar luminosity profile (assuming a 
spatially uniform stellar mass-to-light ratio) and a component 
representing the central point mass. This procedure allows us 
to examine a two-dimensional parameter space of the stellar 
mass-to- light ratio, M/L, an d the central point mass, Mbh- As 
noted bv lBoker et all (119991) . an upper limit to Mbh can be de- 
termined by modeling the cluster in the limit that the stellar 
M/L approaches zero. 

Since the cluster's mass and radius imply a half-mass re- 
laxation time of order ~ 5 Gyr (lMerrittll2008.) . and the cluster 
probably contains significant population components younger 
than this, the cluster may not be relaxed and it could have an 
anisotropic velocity dispersion. When scanning the parameter 
space, we fix the velocity anisotropy and assume it to be con- 
stant throughout the cluster To explore the possible effects of 
anisotropy, we calculate models for three fixed values of the /3 
parameter, jS = -1,0, and 0.5. The values (3 = -I and 0.5 cor- 
respond to ratios of tangential-to-radial velocity dispersions 
of \/2 : 1 and 1 : V2, respectively. 

In Figure[8] we show the resulting two-dimensional param- 
eter space of Mbh vs. stellar mass-to-light ratio. For the spe- 
cific case of isotropic velocity dispersion, the shaded band 
in Figure [8] shows the projected contours of line-of-sight ve- 
locity dispersions corresponding to the la error range on our 
measurement of the cluster velocity dispersion of ct = 43 ± 3 
km s"' . We see from Figure |8] that, in the limit that the stellar 
mass-to-light ratio approaches zero, we obtain a strict upper 
limit to the mass of the black hole, 3 x 10'' Mq. Since we 
are averaging the velocity dispersion over the entire surface 
area of the cluster as in equation |4l this result turns out to 
be nearly independent of velocity anisotropy. We also find 
that, in the limit that Mbh ^ 0, the stellar mass-to-light ra- 
tio is 1 .4 ± 0.2 (M0/L0 y), again independent of the velocity 
anisotropy. 

In Figure |9] we show the two-dimensional parameter space 
of black hole mass and integrated stellar mass within 17 pc. 
In the limit thatM^n 0, the stellar mass within 17 pc is con- 
strained to be within (0.95-1. 2)x 10^ M©. This derived total 
stellar mass is on the larger end o f the distribution of those 
measured bv lWalcher et alj (l2005h . due to the relatively high 
stellar velocity dispersion measured for NGC 3621. We also 
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Fig. 8. — Dynamical modeling results. The curves show the locus of 
models yielding a line-of-sight velocity dispersion of cr = 43 km s"' in the 
stellar {M/L)y vs. Mbh plane. The dotted curves are for radially anisotropic 
models, the solid curves are for the isotropic models, and the dashed curves 
are for the tangentially isotropic models. The shaded band represents the Icr 
uncertainty range on ct* of 40-46 km s"' for the isotropic model. 

note that this result implies a value of / « 6 for the numerical 
factor in the virial theorem mass estimate (Equation[T]i. 

We also computed a set of models in which the outermost 
Gaussian component from the WFPC2 GALFIT decomposi- 
tion (with FWHM = 3"1 or 99.2 pc) was assumed to belong 
to the cluster rather than to the galaxy disk, to assess the im- 
pact that this choice has on the model results. The difference 
proved to be insignificant, because the width of this compo- 
nent is so large that only a small fraction of its mass lies within 
our fiducial 17 pc radius for the dynamical mass calculation. 
In the limit of Mbh = 0, the dynamical mass for the cluster 
changed by < 1 % when the outermost Gaussian component 
was included as part of the cluster. Similarly, the upper limit 
to Mbh is unchanged if this component is added to the cluster 
model, and the eiTor contours in Figures |8] and |9] are essen- 
tially unaffected by this change. 

It bears repeating that in our mass estimates we assume that 
there is no rotational component contributing to the measured 
velocity dispersion. As argued above, this is a plausible as- 
sumption, although it is at this time not possible to determine 
what fraction of the observed velocity dispersion comes from 
ordered rotational motion. Any rotational motion will affect 
the above mass estimates, though determining by how much 
requires a full rotational model of the cluster, which depends 
not only on the rotational speed but also the inclination and 
eccentricity of the cluster If the cluster were rapidly rotating 
and oriented with its rotation axis close to our line of sight, 
then our method would underestimate the true mass of the 
cluster and underestimate the upper limit to Mbh- With the 
data currently available, it is not possible to resolve the in- 
ternal kinematics of the cluster or to constrain the cluster's 
rotation, or to measure the rotation speed or dispersion of the 
inner disk in the neighborhood of the cluster. However, it is 
unlikely that the cluster would be rotating in a nearly face-on 
orientation, s ince the host galaxy as a whole has an inclination 
of / = 65.6° (iPaturel et al.ll2003l) . The high-resolution kine- 
matic st udy of the nuclear c luster in the edge-on galaxy NGC 
4244 bv lSeth et all (l2008bl) demonstrates that rotation can be 
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Fig. 9. — Dynamical modeling results. The curves show the locus of 
models yielding a line-of-sight velocity dispersion of cr = 43 km s"' in the 
stellar mass vs. Mbh plane. Similar to Figure [8] the dotted curves are for 
radially anisotropic models, the solid curves are for the isotropic models, and 
the dashed curves are for the tangentially isotropic models. The shaded band 
represents the Icr uncertainty range on of 40-46 km s"' for the isotropic 
model. 

dynamically important in nuclear clusters, and similar obser- 
vations of the NGC 362 1 cluster will be important for refining 
th e cluste r mass estimate. 

iMerritj ( Il987h has shown that a lower limit to the mass of 
a cluster can be derived from the virial theorem, based on the 
extreme limiting case in which all of the mass is contained at 
the central point of the cluster This minimum mass is given 
by GMniin = 3crJ/(r"'), where l/(r"') is the harmonic mean 
radius of stars in the cluster and ct* is the line-of-sight ve- 
locity dispersion. Averaged over the three-dimensional, de- 
projected light distribution, (r"') = 0.5 pc"', and the resulting 
lower limit to the cluster's dynamical mass is ^ 2.6 x 10^ M©. 
This limit corresponds to our Jeans equation model for a clus- 
ter in which all of the mass is in the central point source, and 
the masses derived by these two methods turn out to be nearly 
identical. We interpret this in the following way. The mini- 
mum possible total mass of the cluster is found for the limit 
in which all of the mass is in a central point, whether calcu- 
lated by the virial theorem or by more detailed Jeans equation 
modeling. As our Jeans equation modeling demonstrates, the 
lower limit to the mass of the cluster overall is the same as the 
upper limit to the mass of a black hole in the cluster. There- 
fore, the quantity 3(T^/G(r"') may be useful as a simple virial 
estimator of the maximum possible mass of a black hole in a 
star cluster This could prove useful as a shortcut for deriv- 
ing upper limits to black hole masses in larger surveys of nu- 
clear clusters, and the limit derived in this way is substantially 
smaller than the total cluster mass derived from Equation[T] 

In principle, the upper limit to Mbh could be lowered fur- 
ther if the stellar mass-to-light ratio could be constrained more 
accurately from multiband photometry and spectroscopy. 
From ^2.31 the photometric data implies a V-band stellar 
mass-to-light ratio in the range (1.4-2.O)MQ/L0y. If we 
conservatively assume that MjLy > 1 then the upper limit to 
Mbh would be ^ 10^ Mq. Since our stellar mass estimates 
from photometry are very uncertain, we keep the extreme case 
of Mbh < 3 X 10^ as our best estimate of the upper limit to the 
black hole mass in this galaxy. 
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5. DISCUSSION 
5.1. AGN Energetics 

The AGN in NGC 3621 is faint and may be heavily ob- 
scured, and it is diffi cult to determine its bolometric lumi- 
nosity. ISatvapal et al. (2007) found a correlation between 
mid-infrared [Ne V] emission-line luminosity and bolomet- 
ric luminosity for a small sample of nearby AGNs with well- 
sampled spectral energy distributions, and used this correla- 
tion to derive an estimate of L[,oi for NGC 3621 based on its 
observed [Ne V] flux. They found Lboi ~ 5 x 10^' ergs s"', 
which is a substantial luminosity for an AGN in a late-type 
spiral. For comparison, N GC 4395 has Lbd ^ 5 x 10"*" ergs 
s"' an d Lboi/i-Edd « 10"^ (iMoranet all 1 19991: iPeterson et all 
l2005h . If the infrared-based estimate of Lboi in NGC 3621 is 
correct, then the AGN must be almost entirely obscured in the 
optical; this would be consistent with its Type 2 classification 
and starli ght-dominated continu um. Combining the value of 
Lboi from ISatvapal et all (l2007l) with our constraints on the 
black hole mass, we have Lboi/LEdd ^ lO""'. 

The [O III] A5007 flux measured from our ESI spectrum is 
(22.4 ± 0.4) X lO""" ergs s"' after correcting for Galactic ex- 
tinction. The uncertainty reflects only the statistical error in 
the extracted spectrum; additional error will arise from the un- 
known difference in slit losses between the gal axy a nd stan- 
dard star observation. As described above in ^3.31 the ESI 
spectrum probably includes only ~ 10% of the total [O III] 
flux from the AGN, so our measurement gives a lower limit to 
th e true emission- l ine lum inosity. 

iHeckman et al.l (l2004t) find that for luminous Seyfert 1 
galaxies, Lboi/L([0 in])f« 3500 with a scatter of approxi- 
mately 0.38 dex. Applying this [O III] bolometric correction 
to the lower limit for the [O III] luminosity of the NGC 3621 
nucleus, we arrive at a lower limit of Lboi > 4 x 10'*" ergs s"' 
based on the optical spectrum. This is an order of magnitude 
lower than the estimated Lboi from the Spitzer spectrum. If the 
ESI spectrum represents 10% of the total nuclear [O III] flux, 
then the full [O III] luminosity would give an estimate of Lboi 
that is essentially consiste nt with the luminos ity derived from 
the [Ne V] emission bv .Satvapal et al.l ( l2007b . Ultimately, the 
best constraints on the AGN energy budget and obscuration 
could be obtained from new X-ray observations, but there are 
none available at present. Perhaps most important, detection 
of a compact, hard X-ray source in the center of NGC 3621 
would provide confirmation that it is genuinely an AGN. 

5.2. Black Holes in Late-Type Spirals 

Our primary result is that NGC 3621, which is a bulgeless 
or nearly bulgeless disk galaxy, contains a nuclear star clus- 
ter in which the maximum possible mass of a central black 
hole is 3 X 10^ Mq. NGC 3621 is one of only a few ex- 
amples of a bulgeless disk galaxy containing both a nuclear 
star cluster and a spectroscopically detected AGN. Under the 
standard assumption that the nuclear activity is driven by gas 
accretion onto a black hole, our results provide new evidence 
that some bulgeless spiral galaxies do contain low-mass cen- 
tral black holes. It is important to acknowledge that the de- 
tection of nuclear activity in NGC 3621 remains a good deal 
less secure than, for example, the well-studied Seyfert 1 nu- 
cleus in the late-type spiral galaxy NG C 4395. Still, the high- 
ionization [Ne V] emission found by ISatvapal et al.l (l2007h 
and the Seyfert-type optical spectrum both provide a reason- 
able basis for classifying this object as an AGN. In this case, 
the low luminosity of the AGN is actually somewhat advan- 



tageous: with no bright nuclear point source it is possible to 
obtain accurate measurements of both the stellar velocity dis- 
persion and the structural properties of the nuclear star cluster 
Such measurements would be far more difficult for a galaxy 
like NGC 4395, which contains a bright, unobscured AGN 
w ithin its nuclear clu ster. 

iRossa et aTl (l2006b found an anticorrelation between mor- 
phological T-type and nuclear cluster mass, and a very weak 
correlation (with very large scatter) between total host galaxy 
B-band luminosity and nuclear cluster mass. Our dynamical 
models imply a stellar mass of ^ 10^ Mq for the NGC 3621 
nuclear cluster if the black hole mass is small compared with 
the cluster's total m ass. NGC 3621 h as T = 6.9 ±0.4 and 
\og{LB/LQ^B) = 9.9 jPafiirel et al.ll2003l) . Wifli these values, 
and assuming a stellar mass of 10' for the NGC 3621 nu- 
cle ar cluster, it i s a fa irly typical object in comparison with 
the iRossa et aTl (l2006l) sample, falling within the scatter of 
both the cluster mass vs. T-type and cluster mass vs. host lu- 
minosity relationships. The only property of this cluster that 
appears unusual relative to previously studied nuclear clusters 
in late-type spirals is its larger velocity dispersion. 

The location of NGC 362 1 in the Mbh - o" relation is shown 
in Figure \TU\ In this case, we use the sample compiled by 
iTremaine et al.l ([2002) for comparison, and include measure- 
ments or constraints on Mbh and a^, for other late-type spi- 
ral galaxies (NGC 4395, IC 34 2, M33, and NGC 1042), for 
the dwarf elUptical NGC 205 dValluri et all lIOOSi). a nd for 
the globular cluste rs Gl and uj Cen dOebhardt et al.l l2005t 
iNoyola et al.l2008h . For the late-type spirals, cr^ is taken to be 
the integrated line-of-sight velocity dispersion of the nuclear 
star cluster The upper limit for NGC 362 1 falls well above the 
Mbh - o" relation extrapolated to low a. Some studies of the 
low-mass end of the Mbh - c relation have suggested a possi- 
ble flattening in slope at the l ow-mass end dBarth et al.ll200l 
iGreene & Hol l2006b;'Wvithe'2006). and it IS important to test 
these claims with new dynamical measurements of black hole 
masses in low-dispersion galaxies, but our results do not pro- 
vide any strong new constraints on the Mbh - o" slope. 

The upper limit to Mbh in NGC 1042 also lies considerably 
above the Mbh - o" relation, but we note that this upper limit 
(from Shields et al. 2008) is simply the total dynamical mass 
of the cluster from Walcher et al. (2005, 2006), which is dom- 
inated by the stars in the cluster. The upper limit could be 
considerably tightened by calculating dynamical models with 
a central point mass in addition to an extended stellar com- 
ponent, as described above. Similarly, it would be interest- 
ing to compute upper limits based on this technique for all of 
the o bjects having s tellar velocity dispersion measurements 
from lWalcher et all ( l2005h . even for nuclear clusters that do 
not contain AGNs, since the more stringent upper limits that 
could be derived from this analysis could provide important 
new constraints on black hole demographics in late-type spi- 
rals. 

With just a few nearby examples of very late-type spirals 
that are known to host AGNs, and just one object (M33) in 
which a massive black hole can be essentially ruled out, we 
have only a very incomplete picture of black hole demograph- 
ics in late-type spirals. An important motivation for improv- 
ing the black hole census in low-mass galaxies with small ve- 
locity dispersions is that the Mbh - o" relation and black hole 
occupation fraction together encode key information about the 
mass scale of black hole s eeds and the e fficiency of seed for- 
mation, as discussed by iVolonteri et al.l (i2008.) . Scenarios in 
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Fig. 10. — The Mbh ~f relation including nuclear clusters in late-type spi- 
ral gal axies. Black circles are objects from the compilation of Tremaine et al.l 
(2002). Filled squares represent nucle ar clusters co ntaining AGNs, includitig 
NGC 3621 (this work), NGC 4395 I Filippenko & Ho 2003; PetersoneLafl 
1200^ . and NGC 1042 ( Walcher et al. 2005). Triangles represent nuclear 
clusters in late-type spirals without AGNs: IC 342 (Boker et al. M33 
{Gebhardt et al. 2001). Globular clusters w ith stellar-dynamical detections of 
black holes are Gl (Gebhardt et alj|2005l) and Cen (Noyola et al. 200^ 
The u pper limit to the black hole mass in NGC 205 is from Valluri et al] 
120051) . Lines denote fits to the Mbh ~ re lation fromlTremaine et al. (2003) 
(solid line) and lMerritt & Ferrares j <200H) (dashed fine). 

which black holes begin as massive (~ 10^ ^o) seeds formed 
by collapse of low angular momentum gas at high redshift 
dKoushiappas et alj 120041 iLodato & Natarajanl l2006h gener- 
ally lead to low black hole occupation fractions in galaxies 
with small velocity dispersions. On the other hand, if black 
hole seeds are formed from the remnants of Population 111 
stars and have masses o f ~ 100 Mp,, then the merger tree 
calculations described by IVolonteri et alJ (l2008h predict high 
black hole occupation fractions in low-mass galaxies, but a 
steep slope for the lower end of the Mbh - c relation, with 
some galaxies hosting very tiny (< 10"^ Mq) black holes. 
While these models still do not incorporate all of the rele- 
vant physical processes involved in black hole growth (such 
as gravity-wave recoil kicks), it is becoming clear that a full 
census of black holes in dwarf and late-type galaxies can 
be a fundamentally important diagnostic of the mass scale 
of black hole seeds at high redshift. Much progress can be 
made thro ugh a combination of multiwavelength AGN sur- 
vevs ( e.g.. iGreene & Hoi 120071: iGallo et al.ll20()8t ISeth et alJ 
I2008al; ISatvapal et al.ll2008l; iDesroches & Holl2008l) and im- 
proved stellar-dynamical searches for low-mass black holes, 
although a full census of the black hole content in low-mass 
galaxies is very far from being feasible at present. 

Another possibility is that black holes in galaxies such as 
NGC 3621 might have formed at a later evolutionary stage, 
within the nuclear clusters themselves. The high densities 
of nuclear star clusters, as well as the lik elihood of ongo- 
ing g a s inflow froni the h ost galaxy disk tSchinnerer et al.l 
120031: iMilosavljevid 12004) . offer suggestions that nuclear 
clusters could be favorable environments for the forma- 
tion and early growth of massive black holes. In suffi- 
ciently dense young star clusters with short dynamical friction 
timescales, intermediate-mass black holes could be formed 
by runaway collisions of m assive stars in the cluster core 
dPortegies Zwart et al.ll2004l) . This process is, however, un- 



likely to operate within typical nuclear star clusters in disk 
galaxies. As noted by IWalcher et al.l (l2005l) . the estimated 
dynamical friction timescale for a 100 M© star to sink to 
the center of a nuclear cluster is longer than the star's life- 
time, and as a resu l t th e runaway process described by 
iPortegies Zwart et alJ (1200 4) would no t occur. The same is 
true of NGC 3621; from Equation 1 of IPortegies Zwart et alJ 
(|2004|) we estimate a dynamical friction timescale of ~ 50 
Myr for a 100 Mq star to sink to the cluster center, which is 
much longer than the lifetime of such a massive star Still, 
more extreme conditions might occur during the early his- 
tory of nuclear cluster formation, and the question of whether 
young nuclear star clusters might be favorable sites for black 
hole formation deserves further theoretical consideration. 

6. SUMMARY AND CONCLUSIONS 

Our main conclusions are as follows. 

1. NGC 3621 has a Seyfert 2 optical spectrum, consistent 
with the AG N identification based on its mid-infrared 
spectrum by ISatvapal et alJ (l2007l) . This adds further 
support for the hypothesis that a black hole is present. 

2. For the nuclear star cluster, we find an effective radius 
of = 4.1 pc, a line-of-sight stellar velocity dispersion 
of CT,t = 43 ± 3 km s"', and an absolute magnitude of 
My = -12.3 mag (uncorrected for internal extinction). 
Simple estimates based on the available HST photome- 
try suggest a stellar mass of (1-3) x 10^ Mq. 

3. From two-dimensional image decomposition of 
2MASS images, we find that the galaxy is dominated 
by an exponential disk, with at most a very faint possi- 
ble pseudobulge component being present. We derive 
an upper limit of = -17.6 mag to the luminosity of 
any bulge or pseudobulge component. 

4. From stellar-dynamical modeling of the nuclear clus- 
ter, we find an upper limit to the black hole mass of 
Mbh < 3 X 10^. This result is insensitive to possible 
velocity anisotropy within the cluster The upper limit 
is based on the extreme assumption of M/L = for 
the stellar population. With improved photometric and 
spectroscopic observations it would be possible to bet- 
ter constrain the stellar mass-to-light ratio and reduce 
the upper limit on Mbh- Assuming that the cluster's 
black hole mass is small in comparison with its stel- 
lar mass, the dynamical models give a stellar mass of 

1 X 10' Mq for the cluster 

These results add to the emerging body of evidence that 
some late-type, bulgeless galaxies do contain low-mass cen- 
tral black holes, and that black holes can be found inside some 
nuclear star clusters. It remains the case, however, that there 
are only a few late-type spirals in which a the presence of a 
massive black hole has been either confirmed (based on nu- 
clear activity) or ruled out (from dynamical modeling). 

Further observations of the NGC 3621 nuclear cluster will 
lead to a more complete picture of the properties of this ob- 
ject. The most pressing observational need is for X-ray obser- 
vations, to provide more definitive confirmation of the AGN 
interpretation and a better estimate of the AGN luminosity. 
Improved measurements of the stellar content of the nuclear 
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cluster can be obtained from high-resolution, high-S/N spec- 
tra extending farther to the blue than our ESI data, in addi- 
tion to improved multi-band HST imaging covering the near- 
ultraviolet to near-infrared. Future extremely large ground- 
based telescopes with adaptive optics may make it possible 
to spatially resolve the kinematic substructure of the nuclear 
cluster, and this could lead to dramatically improved con- 
straints on the black hole mass. 
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